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Abstract 

Focus is on two parties with Hilbert spaces of dimension d, i.e. "qu- 
dits" . In the state space of these two possibly entangled qudits an analogue 
to the well known tetrahedron with the four qubit Bell states at the ver- 
tices is presented. The simplex analogue to this magic tetrahedron includes 
mixed states. Each of these states appears to each of the two parties as 
the maximally mixed state. Some studies on these states are performed, 
and special elements of this set are identified. A large number of them is 
included in the chosen simplex which fits exactly into conditions needed 
for teleportation and other applications. Its rich symmetry - related to 
that of a classical phase space - helps to study entanglement, to construct 
witnesses and perform partial transpositions. This simplex has been ex- 
plored in details for d = 3. In this paper the mathematical background 
and extensions to arbitrary dimensions are analysed. 
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1 Introduction 

Entanglement, a non-classical essential feature of quantum theory, has first been 
recognized in 1935 in the connection with paradoxes. Studying it with mathe- 
matical accuracy began 1964, introducing the Bell states. With the new era of 
proposed applications in teleportation, quantum computing and quantum com- 
munication [BW92j,[B93j it became necessary to use the Bell basis of four "magic 
Bell states" describing the qubit. Now, in the process of extending the concepts 
to systems with Hilbert spaces of dimension greater than two, one faces the prac- 
tical task of defining analogous sets of states. The presentation of all possible 
ideal schemes in [ WOl ] reveals a large field of structures that can be chosen of. 

One item that fits into those schemes is presented in this paper. We were 
motivated to study it mainly out of curiosity on the theoretical side. It has a rich 
structure of symmetries which enable deep concrete investigations on the location 
of the border between entangled and separable states. (Compare |VW00j .) It 
would be no surprise when its mathematical beauty will be reflected in practical 
application. 

In this paper we focus on two parties with Hilbert spaces of dimension d, 
i.e. "qudits". For these two possibly entangled qudits we construct an analogue 
to the well known tetrahedron with four mutually orthogonal Bell states at the 
vertices, [HH96J. This magic tetrahedron includes mixed states, inside lies an 
octahedron of separable states. Each of these states appears to each of the two 
parties as the maximally mixed state of its qubit. Considering the duality between 
maps and states, [ZB04] , these states are related to bistochastic maps, |AU82j 
(but only for bipartite systems). We make some remarks on this duality in the 
concluding Section [91 Now, also for the qudits one might prefer states with this 
property. We call it Locally Maximally Mixed. For qubits, any such LMM 
state can be considered as an element of the tetrahedron, [BNT02] , but for d > 3 
the analogue statement is no longer true. So we perform some studies on LMM 
states in Section [2] and identify special elements of this set. In the following 
sections we then recognize a large number of these special states as included in 
our chosen subspace of LMM. 

For the qudit pair we define and study a special simplex ( "generalized tetrahe- 
dron") W. It has d 2 pure states at the vertices, with specified relations between 
them. Its rich symmetry helps to study entanglement, and it fits exactly into 
the conditions needed for teleportation and dense coding as stated in |W01j . We 
explored this magic simplex for d = 3, [BHN06J. In this paper we bring a detailed 
analysis of the mathematical background. This enable us to extend the study to 
higher dimensions. 

Choose some basis in each factor and define a "Bell state", i.e. a maxi- 
mally entangled pure state, in the Hilbert space C ® C d with the vector 




(1) 



.5' 
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On the first factor in the tensorial product we consider actions of the Weyl 
operators defined as 

W k/ \s) = w k ^\ s -e), (2) 

w = e 2ni/d , (3) 
with the identity 

\s-£) = \s-£ + d). (4) 

The actions of the Weyl operators produce mutually orthogonal Bell state vectors 

\n k/ ) = (# M <g>i)|o 0)0 ). (5) 

The set of index pairs (k, £) is a finite discrete classical phase space: I de- 
notes the values for the coordinate in "x-space" , k the values of the "momentum" . 
Remarks on the relation to the physics of the Heisenberg-Weyl quantization we 
have made for d = 3; details on the mathematics follow in Section HI To each 
point in this space is associated the density matrix for the Bell state, the projec- 
tion operator 

Pk,i = \^k,i)(^k,e\- (6) 

The mixtures of these pure states form our object of interest, the magic simplex 

W = { c k/ p k,i I c M > 0, C M = 1 } • (7) 

As a geometrical object W is located in a hyperplane of the <i 2 -dimensional 
Euclidean space {A = ^2 dk/P k ,i | a k ,e £ ^} equipped with a distance relation 
y^Tr(A — B) 2 . Specifying the origin A = 0, it is also equipped with the Hilbert- 
Schmidt norm v^4^, and the inner product Tr(AB) = a k ,eb k ,e- All this is 
imbedded in the d A dimensional Hilbert Schmidt space of hermitian d 2 x d 2 ma- 
trices. We use this Euclidean geometry for ease of calculations. 

The main goal in this paper is the exploration of the borders of SEP, i.e. 
the set of separable states. We find that the structure of the subset SEPflW, 
the analogue to the octahedron of bipartite qubits, is not quite simple. It is 
not a polytope; but a rather detailed study is enabled by the rich symmetry of 
the simplex W. Using part of it, we determine easily first two polytopes giving 
an inner and an outer fence to the border of SEP. These results, among others, 
appear in Sections [3] and HI 

Symmetry is then studied in detail in Section [5j It simplifies also performing 
the partial transpositions of the states in W. This is discussed in Section [71 So 
we get a closer approximation to SEPflW by studies on PPT, that is the set of 
density matrices remaining positive after partial transposition. Here we refer to 
the Peres criterion [P96J which implies that SEP is a subset of PPT. Furthermore 
the partial transposition maps PPTflW into PPTflW, where VV is another convex 
subset of LMM, also defined in Section [7J The partial transposition maps SEP 
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onto itself, so the cases of bound entanglement detected in W are also cases for 
bound entanglement in W. 

Last, but not least, the symmetry of W can be exploited as the symmetry of 
the set of witnesses [TOO] needed there. Here, in Section EJ we use the mathemat- 
ics of convex cones and their duals. It helps to determine exactly the borderlines 
of SEP. This has been done in |BHN06j for d = 3. Extensions to studies for 
higher dimensions will follow (work in progress). 

This study follows two aims. In the main task of investigations the special 
simplex is constructed and its symmetries are stated. Using these symmetries, 
some details in the structure concerning entanglement are explored. In following 
the second trail we check not to have overlooked anything: the symmetry group is 
maximal, the polytopes are optimal. The proofs of having "best possible" results 
afford some mathematical subtlety. We present these subtle investigations in the 
extra Section [HI 

Various mathematical branches are used: theories of numbers, groups, convex 
sets, matrices and Hilbert spaces. But only some basic facts are needed, to be 
found in any introduction or encyclopedia, as [Sch86j . [V64] . [A42j . |W06j . One 
side effect, which unfortunately makes some pain, is the frequent switching of 
the mathematical points of view. Being too strict on the reference to the context 
would make notations cumbersome and difficult to follow. We try to avoid an 
overburdening with symbols. We refer, for example, to the states with the same 
letters as we use for the density matrices representing them. But we are strict on 
not confusing Hilbert space vectors with states. Big Greek letters denote elements 
of the total Hilbert space, small letters, mostly in the environment | ), are used 
for elements of C d . 

Moreover we simplify the notations, omitting the sign for the tensor product 
concerning the two parties, and write e.g. \(p, s) for an element of C d ®C d instead 
of | ip) <S> \s). Our Unitary Operators U, V, W, that occur in this work acting on 
the global Hilbert space are all of product form. They act locally as U, V, W 
on the first factor, and as U, V, W on the second factor. U\<p,ip) = \U<p,Uip). 
Mostly U = 1, with exceptions in Section [51 

2 LMM states 

Elements of LMM are the states p on the Hilbert space H A <8> T~L B = C d ® C d 
which appear locally to each of the single parties as maximally mixed. The 
partial trace in one factor gives the maximally mixed state u on the other side. 

p A :=Tr B p = u A := ~1 A , p B : = Tr A p = u B := ~1 B (8) 



the community working with operator algebras, such a maximally mixed state is known 
as tracial state. 
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We identify special types of LMM states: Its pure states, isotropic states, 
Werner states and maximally exposed elements of SEPflLMM. 

The pure LMM states for qubits are known as the "Bell states" . We extend 
this naming to each one of the pure LMM states of qudits. The single Bell states 
for fixed d are all unitarily equivalent, involving local unitary transformations: 
Consider the pair Q, $ of Bell state vectors. For Schmidt decomposition, we 
choose a preferred basis {\s)} in Hb- Then there are two different bases \ip s ) and 
\cp s ) in Ha, such that 

s s 

is mapped to |<&) by extension of the local unitary operator 

u = J2\^)(H- (io) 

s 

Mixtures of a Bell state Q with the global maximally mixed state 

« = ^i (id 

define the isotropic states (1 — ct)u + Again all the isotropic states 

with the same a are unitarily equivalent. 

Other special LMM states are the lines of Werner states, related to the lines 
of isotropic states by PT, that is Partial Transposition. See [VW00J for the 
appropriate ranges of the parameters a and other details. They are also all 
equivalent. We need no special check of their belonging to LMM. There is the 
general fact: 

1 LEMMA. The LMM property (GJ) is preserved under partial transposition. 

Proof. The equation (jSJ) is equivalent to the statement that for each <p e C d and 
each basis {ip t } 

^2((p,A\p\<P,^t) = \W\\ 2 /d- (12) 
t 

We write this as TipQ = \\ip\\ 2 /d with Q := J2t \<Pi ■ The PT operator 

in the Hilbert Schmidt space is symmetric, i.e. Tr(PT(p) • Q) = Tr(p • PT(Q)). 
Calculation of PT is done in the preferred basis |s); the expansion \ip) = J2 S <p s \s) 
gives 

PT(Q)=Y,<P.tiPTMt)(r,1>t\) 

s,r,t 

= ^2 ( Ps^* r \r,tpt)(s,tpt \ = W^t){v^tl (13) 



s.r.t 



with \ip) = J2 r f*\r). The complex conjugation does not change the norm and 
(H2D holds for PT(p). □ 
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The Werner states do not only have special symmetries, they have the prop- 
erty of attaining the minimal possible distance (see |GB 02j) to the maximally 
mixed state to when they are at the border between SEP and the entangled 
states. We conjecture that they are the only LMM states with this property. 
What we can show easily is the 

2 LEMMA. If an LMM state at the border between PPT and non-PPT states 
has the minimal border distance to uj, which is l/d\/d 2 — 1, then it is a Werner 
state. 

Proof. Performing PT on the density matrix p of this state we get a density 
matrix o at the border of PPTflLMM to non-positive matrices. The matrices 
with minimal distance at that border have the form a = (1 — P)/(d 2 — 1), with 
P a projector belonging to a pure state. Pure LMM states are Bell states, so cr 
is isotropic and p = PT(a) is a Werner state. The Euclidean distance squared is 
easily calculated as 

^o-^ = ^-^ 2 = ^ 2 -^ = lH^Ty < 14 > 

□ 

W does not contain Werner states if d > 3, but W does; see Section [3 
A third kind of special LMM-states appears in both W and W: The separable 
states with the largest possible distance between uj and SEPflLMM. 

3 THEOREM. The maximal distance of a a e SEPflLMM to u is y/d - 1/d. 
It is attained if and only if the density matrix a has the form 

s 

where both (p s and ip s are bases for C d . 

Proof. The first condition, a gSEP, is fulfilled iff a can be represented as 

A a |y? Q , ip a )((p a , ip a \ (16) 

Q 

with \ a > , J2 a ^« = 1, and normed vectors ip a , ip a . 

The second condition, a GLMM, implies that V<y9, V?/> with norm one 

(tp,^\a\ip,ip) < ^2((p,if>j\<r\(p,il>j) = (<pWa\<p) = 2' ( 17 ) 

3 

where we considered some basis ipj containing the given ip. Applying (fTTj) to the 
vectors appearing in (|T5|) gives 

Trcr 2 = ^2X a (ip a ,ip a \a\<f a ,ip a ) < ^K^i = K (18) 
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This proves the first statement about the maximal distance, since \\a — u\\ 2 = 
Tra 2 — 1/d 2 . To prove the second statement observe that the inequality (TTTj) 
turns to an equality iff VV>~ L with (?/>- L |-?/>) = the equality ((p, ip^^ip, ip x ) = 
holds. The same is true with the roles of the two sides interchanged, that is \/<p 
with (<f \<p) = one has ((p x , iplcr^- 1 -, ip) = 0. So one can start diagonalizing the 
matrix a. One begins with one pair of vectors appearing in ffTTj]) . say a = 0. 

1 . d 1 

a = -\ipo,ipo){¥o,ipo\ H ~(j~ (Jd - 1 - ( 19 ) 

The matrix o d -\ is a normalized density matrix in the LMMflSEP with lower 
dimension, acting on C d_1 ® C d_1 . This can easily be seen by (d — l)TiBO'd^i = 
dTi B a - |y?o)(<A)| = 1 A - |y?o)(y?o| = (d - 1)1^-1 and ((p- L ,ip- L \a d - 1 \ip ± ,ip ± ) = 
djid — l)(v 2 " L , ^Wl^i i J± )- Now one may proceed inductively, expanding c^-i 
in the form (TT6]) - generally with new vectors - diagonalizing (d — l)o~d-i = 
V'll + (d — 2)a d -2, and so on. □ 

Each one of these maximally exposed SEPflLMM states is in unique corre- 
spondence to a pair of bases in the Hilbert spaces of the parties and a one two 
one mapping between them. Each one can be represented as a mixture of d Bell 
states appearing in some W. For example one may use the bases characterizing a 
as stated in Theorem [31 to construct a simplex W: Put them into the definition 
flU), |O 0j0 ) = ■j%Yj a \ l P8,i>s), construct the P kfi , and represent a = Y.k P k,o/d. 
But this representation is not unique. More about this is presented in Section HI 
PT maps this set of maximally exposed states onto itself, so these states appear 
in W also. 



3 Subsets of LMM 

Let us proceede and look at subspaces of LMM. Most important are the Bell states 
appearing in the chosen subspace. Any set of d 2 mutually orthogonal Bell states 
P a - orthogonality of the Hilbert space vectors (fiojf^) = 8 a> p is in this case 
equivalent to the orthogonality of the density matrices in the Euclidean space 
TrP Q P/3 = 5 aj f3 - span a maximal simplex. Each Bell state comes with an optimal 
witness, a hyperplane B a defined as B a := {p : Trp(P a — t/d) = 0}. These 
d 2 hyperplanes B a , together with the d 2 hyperplanes A a := {p : TrpP a = 0} 
containing the faces of the simplex, define an enclosure polytope. Outside 
of it are only entangled states. The projectors P a generate a maximal abelian 
subalgebra of operators acting on TC. So these conditions alone bring already some 
insight, but they still allow for many different choices of an LMM subspace. They 
are not all equivalent. The geometric symmetry of the enclosure polytope, the 
same symmetry as that of the simplex, is deceptive: SEP must be inside, but the 
relations of its detailed geometry to the set of pure states in the chosen subspace 
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depends on their algebraic relations. The single Bell states are equivalent, but 
already pairs of orthogonal Bell states fall into different classes of pairs if d > 4. 
Enter spectral theory: Each class is characterized by the spectrum of the local 
unitary operators identified in ( 1101) connecting the pair. Orthogonality of the 
Bell state vectors implies Till = 0. Being interested only in the intertwining 
relation 

u\o)(n\rf = |$}($| (20) 

we are free to choose a phase factor for U such that one of its eigenvalues is equal 
to 1. The condition Till = specifies the rest of the spectrum only for qubits 
and qutrits. For d > 4 there are various possibilities, defining different classes 
of equivalent pairs: Unitary or antiunitary local mappings of one pair onto the 
other can be applied to the interwiners U. So their spectra are either unchanged 
or complex conjugated and rotated. This characterizes the classes. 

To choose special sets of Bell states an extra criterion which a theoretician 
likes to pose is that the intertwining operators form a unitary group, allowing 
for multiplication of any two of them. This gives a strong restriction on their 
spectra. Enter number theory: 

4 THEOREM. If {U n } is a group of intertwiners between mutually orthogonal 
Bell states, then, with an appropriately chosen overall phase factor, U has eigen- 
values e 2lTim / b j where < m < b — 1, and b is either a divisor of d or equal to 
d. Considering intertwiners acting as U = U ® 1 on Ti^ only, the multiplicity of 
each eigenvalue of U is d/b. 

Proof. The Euclidean space of density matrices has finite dimension, the set of 
orthogonal projectors onto U n \Q) is finite, less than d 2 , and there exists some 
smallest natural number b, such that U b = 1 ■ phasef actor. We choose the phase- 
factor for U in such a way that we have U — 1. In the following we consider 
U acting on T~Ca only. Since Trf/ 6 = d Trf/ 6 , the orthogonality of the Bell states 
implies, as stated before equ. (T2"U]) 

TrU n = d5 nfi for < n < b - 1. (21) 

U b = 1 implies that the eigenvalues of U are elements of { e 27rjm / fe j < m < 6—1} 
Denote the multiplicities as f{m). Then the equation ( 1211) can be read as a 
formula for the Fourier transform of f(m). The inverse transform gives 



1 6-1 d 
f(m)= J2e 2nimn/b TrU n = j-. (22) 



b^ b 

n=0 

This number has to be an integer. □ 

5 COROLLARY. Any group of unitary intertwiners between mutually orthog- 
onal Bell states contains finite cyclic subgroups. Each one is of some order b, 
where either b = d, or b is a divisor of d. 
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It follows that there are not many different possibilities for structures of such 
groups. Our choice is possible for all d, whether prime or not. 

4 Groups and the classical phase space for the 
magic simplex 

Letters of the set {j, . . . t} denote numbers 0,1 ... d — 1. They are considered 
as elements of Z^ := Z/cfZ. Calculations with them are to be understood as 
"modulo d". 

Intertwiners are the Weyl operators W^/ — Wk,i <8> 1 presented in Section [T] 

WkjPp,qW k j = P p+ k,q+£ ■ (23) 

The Weyl operators obey the Weyl relations 

W jtt W k , m = w ke W w+m , (24) 

Ke = W k~J = ™ M W- k ,-e, (25) 
W 0fi = 1. (26) 

They form the Heisenberg-Weyl group W. More precise: W is a finite discrete 
subgroup of the doubly infinite continuous Heisenberg group; compare [W06J. 
Group elements are w m Wk/. The phase factors {w m l} form an abelian normal- 
ize^ the factor group is W/Z^ = Z^ x Z&. This can be considered in the sense 
originally meant by Weyl, [W31] , as the quantization of classical kinematics. The 
kinematics of the Galilei group is represented in the discrete classical phase space 
as Zd x Zd, generated by the global boost (p, g) i— > (p + 1, q) and the global space 
translation (p, q) \— > (p, q + 1). 

The classical phase space T := {(p, q)} is a lattice on a two-dimensional torus. 
It has a "linear" structure - multiplication by constants and addition is always 
done in the ring Z^ - and it is a symmetric space for the Heisenberg-Weyl group: 
We define the action of W on T by identifying each phase space point (p, q) 
with the projector P pq and use equ. (!23|) . Moreover we identify non- negative 
normalized densities {c Pi9 > 0, Yl c p,g = 1} with the elements ^c p>q P p ^ of W. 
Special use is made of equidistributions over subsets Q C T and the corresponding 
density matrices 

Now the group structure of W gives a first insight into the structure of SEPflW. 
Each cyclic subgroup {W^} acting on a point (p, q) of T generates a line 
{(p + nk, q + n£)}. □ These lines, for d prime, have been identified in [N06J 

2 Warning: These lines do not for each d fulfill the conditions for a "line" in the sense of 
affine geometry. See also |B04j . 
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as corresponding to separable states. If there are non-cyclic abelian subgroups 
- which may be the case if d is not prime - they generate sublattices, each one 
with at most two independent basis vectors. 

6 PROPOSITION. Each line or sublattice with d points is generated by an 
abelian subgroup o/W and corresponds to a maximally exposed state in SEPdW. 

Proof. Consider a sublattice Q with d points. A lattice in 2 dimensions can be 
represented with 2 basis vectors. So we can represent 

Q = {(b + j/i + kv,q + £fi + mis), < fi <b - 1,0 < u < c-l,b - c = d} (28) 

We include the cases b = d, c = 1, representing lines. For the matrices the 
representation is 

PQ ^^EE^I^fel^" (29) 

where U = e n Wj/, V = e l5 Wk, m , with the phase factors chosen, if necessary, such 
that U b = V c = 1. Each one of the smaller exponents gives other elements of 
W; so 

TtaWV = ■ d (30) 

For the sublattice the Weyl relations < HM imply 

U ■ V ■ C/ f ■ = w M - jm l, (31) 

and the exponent k£—jm is the oriented area of a unit cell of Q. The union of all d 
cells spans all of T, which has area d 2 , once or several times; so (k£—jm)-d = z-d 2 , 
with some z6Z. It follows that k£ — jm = 0, the r.h.s. of (!3~TTU s 1, so U and V 
commute. This allows for a common spectral decomposition 

d/b-l d/c-1 

V = E E /M)e 2 " s/ WKtl, 

s=0 t=0 
d/b-l d/c-1 

V = E E f{s,t)e 2 ™ t/c Ws,t)M. (32) 

s=0 t=0 

We get the multiplicity function /, in a way analogous to the proof of Theorem 
HI Here we use the Fourier transform in two variables, and equ. (l30|) : 

f(s, t) = — e2m W+rtrtTrf/rf" = 1. (33) 

The diagonalizing basis (p s ^ in 1~La is now used for a Schmidt decomposition of 
the Bell state vector. 
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With ip s j as the appropriate basis in Hb- Inserting (132]) to (|34|) into (|29|) . the 
summation over the phase factors brings some S factors, reducing the summations. 
The result is 



This expression for pq is exactly as it is used for the matrices in Theorem □ 

SEP is a convex set and the separable states pq identified in Proposition [6] 
can be considered as the extreme points of a kernel polytope which is a subset 
of SEPflLMM. These pq appear also as extreme points of the enclosure polytope, 
but do not cover all of them if d > 3. For the vertices of the enclosure polytope 
the set Q can be any subset with d elements of T; for the kernel polytope this 
set Q has to be a line or a sublattice. That all the other sets Q correspond in 
fact to entangled states is proven in the Theorem [14] in Section [HJ 

7 THEOREM. The number of lines and sublattices with d points in T is 



where the sum runs over all the b which are proper divisors of d. 

Proof. The number in square brackets must be the number of lines and lattices 
Q of order d, each containing the point (0, 0) G T. All the others can be found 
by translations; and doing all d 2 translations gives each line and each lattice of 
order d in d-fold multiplicity 

We give a list of these Q containing (0, 0): 

a) {0,0)| 0<s<d-l}, one line 

b) {{k ■ s, s)\ < s < d — 1}, d lines, one for each k e [0, . . . d — 1] 

c) {(n ■ b + v ■ b, v ■ d/b)\ < p, < d/b - 1, < v < b - 1}, b sublattices, 
one for each v e [0, b — 1], with b a proper divisor of d. 



We remark that some of the sublattices listed in c) can also be considered 
as lines. But not all of them, if d is not a simple product of prime numbers but 
contains also squares or higher powers of them. One example is d = 4, b = 2, 




(35) 



s.t 




□ 



v = 0. 
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5 Symmetries of W 

We are looking for symmetries compatible with the entanglement, just to make 
the investigations simpler. We do not pose detailed restrictions, no measure for 
entanglement is needed. Just the following, physically motivated characterization 
is sufficient: 

8 DEFINITION. A mapping L : W —> W is E-compatible (i.e. compatible 
with entanglement), iff 

a) Bell states are mapped to Bell states, 

b) L is mixture preserving 

L{ap + (1 — a)a) = aL(p) + (1 — a)L(a), 



c) SEPPiW is mapped onto itself. 

Now if we have a local unitary transformation p — > UpU* , the separability is 
preserved. Also the conditions a) and b) are fulfilled. So we know already about 
a subgroup of symmetry transformations: The translations of phase space; they 
are implemented as local unitary transformations in the Heisenberg Weyl group 
W. 

For general L the fulfilling of b) implies the possibility to linearly extend L. It 
gives then, due to a) , an invertible norm preserving linear map of the Euclidean 
space, spanned by the Pk/, onto itself. It effects a permutation of the Euclidean 
basis elements Pk/, hence a map T — > T. Vice versa, any permutation of this kind 
extends via mixing preserving to a map W — > W. Now, after any permutation 
of T a certain translation can bring the origin (0, 0) back to its place. So each of 
the symmetry operations can be formed as a product of a phase space translation 
with a certain point transformation M, which leaves the point (0, 0) at its place. 

The tool box of these point transformations of phase space contains: 

The "horizontal" shear of phase space , H : ( ^ ) i— > ( | . Its 

; \ q J \ p + q J 

powers form a cyclic subgroup. The elements are represented by the matri- 
ces ^ ^ ^ ^ . (In discrete classical mechanics this is a free time evolution.) 

The "vertical" shear of phase space ,V: (?) ^ (^?^)' ^ 6 
elements of the generated cyclic subgroup are represented by the matrices 
( 1 ) ' ^ m8 ^ ^ e cons idered as a local boost.) 
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A quarter rotation of phase space , 1Z : ^^Ji-^^^J.Itis rep- 
resented by the matrix ^ ^ y ' ^ = — ^ s a P om t reflection. 

Squeezing , a scale transformation ^^^h^^ r ^ j , for r • s = d + 1, 
possible for each r relative prime to d. 

Reflections : 

• Inversion of momentum, S : p i— > — p 

• Space reflection, g i— » —q 

• Diagonal reflection, p i— » g, g i— > p 

9 PROPOSITION. Consider a linear mapping T — > T, defined as the appli- 
cation of a 2 x 2 matrix M with elements G Z^ ; and wi/i det M = ±1. 5y 
extending it to a mapping W — > W, is E-compatible. 

These matrices form the extended symplectic group Sp(2, Z<$). 

Proof. Addition and multiplication of the matrix elements is according to the 
rules of the ring Z^. This defines the matrix multiplication. The unit matrix has 
det 1 = 1 and is also in this set. Inverting a general element M is achieved with 
the mapping 

f k m \ „.i . / n — m \ 

M =(< „) " M =± (-< k j' (36) 

with the sign equal to the sign of det M. So the matrices which are either sym- 
plectic, det M — 1, or mirror symplectic, det M — — 1, form a group. We establish 
now the three transformations V, 5 as generating elements: 

First, they generate 7i = Hr x VR. and all the powers V*, Ti 1 . The space reflec- 
tion is IZ^SIZ, diagonal reflection 1ZS = — 7£ _1 <S, and squeezing is 1ZV~ s 7l r V s . 
Multiplication by S maps symplectic to mirror symplectic matrices and vice versa. 
Consider a general symplectic M. The condition det M = 1 can only then be true 
if k and m are relative prime. Also k and £ have no common divisor. (Letters 
denoting the matrix elements are placed as in (|36l) .) So there exist a t and an s, 
such that m + k t = and £ + k s = 0. One calculates 

\ £ n + t£ J \ n + tt 

This matrix performs a squeezing. It can be represented as stated above. Multi- 
plication by lir s from the left and 7i~ t from the right gives us back the matrix 
M. 
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To realize the E-compatibility we proceed as we did for d = 3. A general group 
element M can be considered as a product of the three generating elements. For 
these we present the operators C,U-]z,Uy and use their products as Um- The 
E-compatibility of the generating elements infers so the E-compatibility of M. 

Now we construct for each generating group element M an operator U m®U^, 
either unitary or anti-unitary. Its local action in the factor C d on the left hand 
side as U M transforms the Weyl operators, unaffected by U M , as 

U M W k/ U^ = e^ M ' k ^W k/// , (37) 

when M maps (k, £) h- > (k', £'). Some phase factors e lv may appear. Then we use 
operators U acting in C d on the right hand side. They are uniquely defined by 
the condition that U ®U~ l leaves the chosen Bell state vector f2 ,o invariant. Its 
matrix elements in our preferred basis are (s\U\t) = (t\U\s). The joint action in 
the space W can now be calculated to give 

P k ,i -> (UM®U M 1 )PkAUM^UM) = U M Pk,iU M 1 =P k ^e. (38) 

Now we look at an implementation of the generating elements as local trans- 
formations of the Hilbert space. The reflection S can be implemented by complex 
conjugation in the preferred basis 

s s 

This is a local anti-unitary operation. It acts onto the Weyl operators as 

CW k/ C = W_ ktl . (39) 

Its anti-linear extension C is complex conjugation in the global Hilbert space, 

CW k/ C = W. k/ , (40) 

mapping SEP onto SEP. So S is E-compatible. 

The other two generators are implemented by local unitaries, so the E-compatibility 
is obvious. The quarter rotation 1Z is implemented as Un = U-r <E> U^ 1 by the 
local Fourier transform: 

t 

It acts onto the Weyl operators as: 

U n W k/ Ul = w- ke W e , k . (41) 

For implementing the vertical shear of phase space, V, one may choose any integer 
v and define 
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For general dimension d we use the ordinary integers s G {0, 1 . . . , d — 1} when 
calculating the exponents. (In Z^ dividing by 2 is well defined for odd d only.) 
For even d the half-integer powers of w have to be chosen consistently for all the 
odd s. This choice, e.g. w M//2 = e l7Tfl ^ d , appears then also in the action onto the 
Weyl operators: 

U v W k4 Ut = v fl*+»*>)fi W w . (42) 

□ 

Now this group of E-compatible point transformations is maximal, other 
transformations do not have the compatibility property. The proof is given in 
Section [HJ 

Some remarks on more subtleties: Note that the implementations of the ma- 
trix group elements M are not unique. There are four different groups of trans- 
formations of vectors and operators involved. Transformations of vectors in C d 
by Umi unitary transformations of the operators as noted in (13"T|) . Then there 
are transformations by Um ® Um °f vectors in the global Hilbert space, and the 
related transformations of the operators as noted in (|38|) . Only the last one gives 
a representation of Sp(2, Z d ), when restricted to the Euclidean space spanned by 
the P k ,i- The others are "quantizations" , involving phase factors. There is more- 
over the possibility to multiply each Um by some Wj jr , and this gives a discrete 
set of different implications. 

The total group E of E-compatible transformations has, by the way, the struc- 
ture of the semidirect product Sp(2, Z^) x W of the extended symplectic group 
and the Heisenberg-Weyl group. That is (see, for example [W06] ) : The extended 
symplectic group gives an automorphism of the normal subgroup W, as noted in 
(J3"7j) . The group product in E is given as 

(U M Wj, m ) * (U L W k4 ) = (U M -L ^ j ' m) W r+k , m , +t ) (43) 
when L maps (j,m) i— > (j',m'). 

6 Witnesses and more symmetry 

Entanglement witnesses have been introduced, |T00j . to detect the entanglement. 
Detection can be either experimentally or theoretically. To prove Theorem [J4] 
in Section [S] we use them in that way, to discern the entangled vertices of the 
enclosure polytope from the separable ones. 

Here we reverse the point of view. "Witnesses" are used to study the loca- 
tion of SEP, the convex set of separable states. We define the set of structural 
witnesses, 

SW := {K = ^ | V<7 G SEP : Tr(aK) > 0} . (44) 

3 Some of the unitary transformations of operators appear as "gates" in Quantum Compu- 
tation, see e.g. [G99] . [GKPOlj . 
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This set forms a convex cone of operators. SWU{0} is the dual convex cone 
to {ap, a > 0, p G SEP} and thus completely characterizes the location of 
SEP. Geometrically, every structural witness defines a hyperplane in the Hilbert- 
Schmidt space of hermitian matrices p, which is a Euclidean space with dimension 
d 4 . The extremal rays of this dual cone are tangential witnesses for density 
matrices p on the surface of SEP. 

TW:= |J TW P (45) 

p£ surface (SEP) 

p G surface (SEP) : TW P := {K G SW | Tr(pK)=0}. (46) 

Being interested in SEP restricted to a linear subspace of states, we may 
restrict the study of witnesses onto a dual subspace. If the set of states is defined 
by invariance under the action of a group Q, the dual subspace is a set of witnesses 
which are also invariant. The details of this argument have been presented for 
d = 3. There was no use of a special dimension and we may take over the results 
from [BHN06| : 

10 THEOREM. Characterizing SEPnW through witnesses is simplified by us- 
ing the following properties: 

• SEPnW is completely characterized by duality, using witnesses of the form 
K = J2k,e K k,tPk,e- 

• Such an operator K is a witness, iff Wtp G C d the operator 

^Kk4W k4 $)$\Wj$ (47) 
k,e 

is not negative. 

• K is a tangential witness in some TW P iff 3\(p,ip) such that 

J2 K k,iMW k M\ 2 ^0 (48) 

k,l 

with \i>) = Y.s( s \^)*\ s )- The state 

P={WA)M\)q (49) 

is a boundary state of SEP and located in the tangential hyperplane, TrpK = 
0. Here Q is the abelian group of unitaries generated by the (1 — P^)/2. 
( )g denotes symmetrizing by the "twirl" operation concerning the group Q . 
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Some states have more symmetry and the tangential witnesses can be found in 
some even smaller set, showing the same symmetry as the state. Sometimes these 
extra symmetries are given as subgroups of the inner symmetries of W which we 
analyzed in Section[5] The simplest example concerns the isotropic witness, which 
is the optimal entanglement witness for a Bell state. The elementary calculation 
may again be performed for general d as it is done for d = 3. Sometimes external 
groups, mapping part of W to other states, have to be used. We give one example. 

11 THEOREM. Concerning the subsection {p = Ylk CkP>k fi} of density matri- 
ces, the search for witnesses can be reduced to {K = J2k K kPk,o + J2e lz Qt}> w ^ 
Qt = J2k p k,i/d, and 7_ £ = j t . 

Proof. As a first step we use the W-symmetry of space reflection i <-> —£. This is 
an invariance of the chosen states. Projecting K = t Kk/Pk,e onto an invariant 
operator by the "twirl" operation with this group Q of two elements gives 

( K )s = K U P W ' K 'k,t = \(. K Ki + (50) 
k,e 

In the second step we use the group QU of local unitaries U ®U diagonal in the 
preferred basis, 

U : \s) i-> e ia{s) \s), U : \t) i-> e~ ia ^\t). 
We use the expansion 

p k,e = ^ fc( *~ r) l* -e,t)(r-t,r\, (51) 

t,r 

and form the second projection by twirl onto operators invariant under this group, 
{(K)g)gu = £ 4,^ E wHt ~ r) (I* - ^ *> < r " ^ r I W ( 52 ) 

The invariant part involves 

(\t-£, t) (r - £, r\ ) ef/ = ( e <WH-«W-«(r-0+«{r)) ^ _ ^ _ ^ r | 
= &,o|t,t)(r,r| + (1 -5«, )5t,r|t-^t)(t-^*| 

Inserting this equation and also (1301) into (1521) gives for the first term J^ fe Kk,oPk,o, 
for the second term (1 - 5 efi heQe with 7^ = ~ ^ki^e + K k) -e). □ 

7 Partial transposition 

PT can be used, referring to the Peres criterion, to prove entanglement. On the 
other hand it maps LMMflPPT onto itself, see Lemma [IJ There is a PT related 
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subset VV of LMM with VVHPPT= WnPPT: It is defined as the linear extension 
of PT(W fl PPT) to the borders of positivity. The dimensions of these related 
subspaces are equal, dim(W) = dim(W) = d 2 — 1. Studies on the structure of 
W are automatically studies on the structure of VV. The two pictures Fig. 2 and 
Fig. 3 presented in [BHN06J for d = 3 can be seen in that way. The region of 
PPT-matrices (not necessarily positive) becomes the region of states, i.e. posi- 
tive matrices, and vice versa. Their intersections are the PPT-states - density 
matrices which are both positive and PPT - in both points of view. The Peres 
criterion, SEPcPPT, implies that also WnSEP= VVflSEP. The cases of bound 
entanglement, [HHH98J, may therefore also be seen in two ways. The regions of 
bound entanglement in W, e.g. those that we found for d = 3, are in one to one 
correspondence to those in VV. 

PT of our simplex VV has nice features, inferring simplification for calculations. 
We use again the expansion (|5"T|) and observe that PT maps 

\t — £, t)(r — £, r\ i— > |m — t, t)(m — r, r\ with m = t + r — £. 

Splitting the global Hilbert space into subspaces according to the quantum num- 
ber m allows for a splitting of partial transposed W-states: 

PT : p = J2 c k/Pk,e >-> 05 m , (53) 

k,i in 

with hermitian d x d matrices B m , 

(s\B m \t) = ~ d Y. c^s+t-mw^ = (t\B m \s)* . (54) 
k 

12 THEOREM. Consider the matrices B m corresponding to some state in W 
according to fiffl). 

• For odd d all the B m are unitarily equivalent. 

• For even d there are two classes of mutually equivalent B m , one for even 
m, the other for odd m. 

• If d is even, there is the relation of matrix elements for every B m 

(s + d/2\B m \t + d/2) = (s\B m \t). (55) 

Proof. For any d observe 

(s\B m _ 2 \t) = ^Y, C k,s+t-m+2W k(s - t) = (S + l\B m \t + 1) 
k 

For d odd one shows (s\B m ^i\t) = (s + (d + l)/2\B m \t + (d + l)/2) by observing 
s + 1 — m + 1 = [s + (d + l)/2] + [t + (d + l)/2] — m in the second index of c. For 
even d, the equivalence s + t = [s + d/2] + [t + d/2] implies (|55|) . □ 
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The last point has the consequence that, if d is even, each B m has the form 
of a block matrix 

£ ^ J ^ C ® 1(2) + D®<r x 2i — ±— © (56) 

with hermitian blocks C m and D m . 1(v) is the v x z/ unit matrix. 
Consider now the abelian algebras 

A{d) := {J2 ak,eP k ,i, a M ^C} = M (d 2 , C) (57) 

emerging as a linear span of the special density matrices. Using Theorem [Tj2] and 
( 156)) the results of mapping by PT are the following subalgebras of M(d 2 , C): 

• If d is odd: PT: A(d) ^ M(d, C) (8 1(d), 

• If d is even: PT : A{d) h+ Af (d/2, C) ® M (4, C) ® l(d/2). 

A consequence is a simplification for checking whether a state in W is PPT or not. 
These states are mapped to linear functionals of PT (A(d)), represented either, 
if d is odd, by hermitian matrices in M(d, C) or, if d is even, by four hermitian 
matrices in M(d/2, C). 

A further consequence is an insight into the structure of the state space W: 

13 THEOREM. The subsetVV of LMM is given by the intersection ofPT(A{d)) 
with the set of density matrices. 

Only for d — 2 it is again a simplex - the reflected tetrahedron. For odd d it is 
the state space consisting of hermitian dxd density matrices - when the tensorial 
factor 1 is neglected. For even d > 4 there are three-dimensional sections with 
the form of a tetrahedron through every point in this d 2 — 1 dimensional convex 
body. In other directions there exist sections of dimension d 2 /A — 1 with the 
structure of the state space with (d/2) x (d/2) density matrices. The space of 
states for M(i>, C) is the convex set of normalized positive v x v matrices. Every 
maximal face is equivalent to the set of normalized [y — 1) x [v — 1) matrices. 
So its faces have dimension v[y — 2) at most. It follows that the surface of W 
is curved in many directions, if d > 3. Part of the border of W is the border of 
PPTflW. This border is therefore also curved in many directions. 

Both local unitary transformations and the global complex conjugation map 
PPT onto itself. So the symmetries established in Section [5] are symmetries of 
PPTflW and of VV too. Also the witnesses for W can be transported to witnesses 
of VV by PT. This follows from the "self-adjointness" of PT as a transformation 
in the Hilbert-Schmidt space: Tr(PT[p]iT) = Tr(pPT[fT|). 

Finally a look at special LMM states in VV. No Bell states are in VV, if d > 3. 
There is a set of Werner states instead, and W includes as many Werner states 
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with some given mixing as W contains Bell states: d 2 of them are extremal with a 
density matrix which is a d(d— l)/2 dimensional projector. The set of maximally 
exposed LMMflSEP states are mapped by PT onto itself. Their number in W is 
thus again N(d), the same as in W, see Theorem U\ and Thm. [THin Section [HI 

8 Optimality 

We follow the second trail which aims at proving not to have overlooked anything. 

14 THEOREM. There is a one to one correspondence between the maximally 
exposed states in SEPDW and the lines or sublattices with d points, generated by 
abelian subgroups ofW. 

Proof. One half of this theorem is proven in the Proposition [61 On the other 
hand, SEPflLMM is inside the enclosure polytope. In the large space of hermitian 
matrices the extremal points of this polytope lie at the intersections of the witness 
hyperplanes B a and the positivity borders Ap, with a G Q C T, (3 G T\Q. 
Restricting the space to the space of normalized matrices gives the condition 
| Q\ = d. This is the condition to get those vertices of the enclosure polytope 
which are inside of W. They all have exactly the same distance to u as the 
maximally exposed separable states. But not all of them are separable; only 
those, where Q is a line or a sublattice. For d prime this has been stated in 
[N06J. For general dimension d we define 



and claim that it is an entanglement witness if e is small and if Q is not a line or 
a sublattice. To prove this claim we have to show that V \(p, ip) the expectation 
value of ([58]) is not negative, 



With P a = Wgn 0j o>(tto,o|Wt one gets 

M\KW^) = IMI 2 IHI 2 -(l+e) \M\w a \si Qfl )\ 2 . (60) 



K =l-(l + e)J2P a = l-(l+e)d-p Q , 



(58) 



(59) 



a 



We insert the definition (CD) of Sl ,o; 



(<^|Wgfi ,o> = ^ d J2&\W a \s)(iP\s) = j= d { V \W a \4>) 



with := Yl s ( s \' l l ; )*\ s ) ■ S° we have to check the non-negativity of 




(61) 
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Since Q is not a sublattice, the Weyl operators which appear in the sum do not 
all commute with each other. That means W a Wp = e n WpW a with e' 7 7^ 1 for 
some pairs of operators, and there is no common eigenvector. For each pair of 
vectors there is at least one a such that | (^W^l^) | < ||y||||'0||- There is only 
a finite number of operators and a compact set of normalized vectors; one has 
equicontinuity and uniform boundedness, 

3e > 0, s.t.V^-0 : ^ \(<p\W a \ij>)\ 2 <d-(l- 2e)|M| 2 |M| 2 . 

So the claim that K defined in (|58|) is a witness for some e is proven: 

(ip^\K\^)>\\ip\\ 2 Uf{s-2s 2 ). (62) 

Since TiKpQ = —e, the state pq is shown to be entangled. □ 

Remark: The procedure connecting the expectations fIBUl) with the formula 
f lBTj) is used also in Section [HI Theorem (TDJ 

In Theorem [TH we have proved that the geometric symmetry of the kernel 
polytope is smaller than that for the enclosure polytope. One implication is: 

15 LEMMA. Every E-compatible point transformation M must be a linear in- 
vertible mapping T — > T. 

Proof. The set of kernel vertices has to be mapped onto itself. This set corre- 
sponds to the set of lines and sublattices with exactly d points in the phase space 
T. Every pair of phase space points lies on one line at least, many pairs on not 
more than one. Since the mappings are one to one, each of these one-line-only 
pairs has to be mapped onto an equivalent one-line-only pair. There are enough 
of them, like [{p, q),(p + k,q + 1)], to imply the linearity: every line is mapped 
onto a line. □ 

We remark that invertibility of M means that the Matrix 



has to exist. This is only then the case if det M is coprime with d, excluding e.g. 
det M = 2 for d = 4, and det M = ±2 or 3, for d = 6. 

Next we show that the geometric symmetry of the kernel polytope is still 
deceptive, if d — 5 or d > 7. 

16 THEOREM. Consider a linear mapping T — > T defined by applying a 2 x 
2 matrix M with elements G Z d . Its extension to a mapping W — » W is E- 
compatible if and only if det M — ±1. 
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Proof. The first part is proven constructively in the proof of Proposition [91 To 
prove the other direction we use duality of convex cones. A linear mapping 
W i— > W is E-compatible iff the dual transformation maps SW to SW and TW 
onto TW. The dual transformation, acting onto witnesses, is given by the dual 
mapping of the set {nk,t} considered as an element of £ 2 {T, R): 

K = } j K kji P k> e, P = ^2 C W P W ^ TrKp = « fej ^ c M , 

The dual mapping of T is therefore M _1 , which is an element of Sp(2,Zd) iff M 
is such a matrix. 

Consider now a line of tangential witnesses 

K(e) = \(e)l + H + eP (63) 

The parameter A is fixed through the conditions on K stated in Theorem [TUJ 
They imply the existence of normed vectors \<f,i(j) such that 

(<f,ip\K\ip, ijj) = min(x, rj\ K \x, v) = 0, (64) 

and therefore 

- A(e) = mm(x,v\{H + eP)\x,v)- (65) 

x,v 

This situation is treated perturbatively. Let \(f(e),ip(e)) be a differentiable curve 
of vectors with \(p(Q),if)(Q)) = the minimizers at e = 0, with normalized 

vectors ip(e) = if + edip + 0(e 2 ), ip(e) = ip + edip + 0(e 2 ). With 

-A*(e) = M£),^)|(tf + £P)|^),V^)> 

one gets 

- ^/x(e) | e=0 = (<^|P|<^) + [<M#» +cc] + [(S^H^ + c.c], (66) 

with the operators and if^ defined as quadratic forms in the local Hilbert 

spaces, 

(x\H^\v) ■= (x, ^|#|»7, V>> , (xl^l 7 ?) := (<P, x\ H \(f, rf) . 

We know from standard perturbation theory that the terms in square brackets 
in (1661) are zero if the ground states of the local and are not degenerate. 
For e = 0we choose 

H = Y1 7fcPfc >° lk = ~{w k + w~ k ). 

k 
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Using the expansion (!5T|) and then ~ ^klkW^ *-* = S S)t ~i + <5 S)t+ i we get 

((p,l/j\H\(p,tlj) = -^2<Pe^e<Pl^t{ S s,t-l + Ss,t+l) (67) 

s,t 

= -2j2mns+i), (68) 

with f(s) := (p s ip s . The minimum of (1681) is attained if all f(s) are real valued and 
positive. Also the tp s and ip s can be chosen as positive. One step in minimizing 
(I6T|) with the condition ||<^|| = ||^|| = 1 can be considered as equivalent to the 
reverse task of keeping the f(s) fixed, and minimizing \\(p\\ ■ This gives 

ip s = ip s and Yl s f( s ) = IMI 2 = !■ Using this as a side condition to minimize 
( l68l) one gets the minimizers: for d > 5 they are f(s) = \(5 s ,t + 8s,t+i) and 
f(s) = \8 Sjt + l(5 Sjt -i + 8 s ,t+i) for any t. Defining ip = v7 and ^ = v 7 / one 
sees the non-degeneracy of the ground states of the local operators and H v . 
Using the ground state vectors v 3 s(^) and ip s (s) gives (i(s) = \(e). Applying (I6B1) 
results therefore in 

-^X{e)\ e=0 =((p,^\P\i P ,^). (69) 
Choose if s = ip s = ^75(^,0 + ^s,i)' consider P = P 0>1 and transformations by 



M 



1 

n 



H(0) is invariant under this transformation, but the perturbative Po,i changes to 
Po, n - The value of ( J69l) is changed, unless n = ±1, proving the non-invariance 
of TW. Using transformations by symplectic matrices, every matrix M can be 
transformed to this diagonal form without changing its determinant. So det M = 
±1 is a necessary condition. □ 



9 Summary and Outlook 

The article extends on previous work [BHN06J for qutrits, but here the results are 
stated in a more mathematical context and are generalized to arbitrary Hilbert 
space dimensions d. We consider the state space of two qudits and analyze a cer- 
tain subset, the simplex (generalized tetrahedron) W as the main object of our 
investigations. It is obtained starting from a certain maximally entangled pure 
state, a Bell type state. By applying on one side the Weyl operators other orthog- 
onal Bell type states are formed, and the set of mixtures, the complex hull, is the 
simplex W. The Weyl operators are related to a discrete classical phase space, 
representing in turn the algebraic relations of the Weyl operators. This analogy 
enables us to describe the local transformations of the quantum state space of 
interest and is very useful for several proofs in this paper: Transformations of 
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W onto itself can be considered as transformations of the discrete classical phase 
space. Thus the symmetries and equivalences can be studied by this means. 

The simplex W is embedded in a d 2 -dimensional Euclidean space equipped 
with a norm (the Hilbert-Schmidt norm) and an inner product. We analyze in 
detail the properties how it is embedded in the whole state space of two qudits 
and discuss symmetries and equivalences inside the simplex W, its facets and 
witnesses. This is obtained via the Weyl group which is a kind of "quantization" 
of classical phase space. 

Then we investigate the question of the geometry of separability. We start 
with the construction of two polytopes, an inner (kernel polytope) and an outer 
(enclosure polytope) fence for separability. They define entanglement witnesses 
but are in general not optimal. The outer fence, the enclosure polytope, has the 
same geometric symmetry as W. Because we are able to construct optimal entan- 
glement witnesses explicitly we obtain in principle the border between separable 
and entangled states, sometimes even in analytic form. 

With our method we find also the set of bound entangled states of the param- 
eter subspace under investigation, by applying the partial transposition on one 
subsystem, which detects entanglement via PPT. The obtained PPT-witnesses 
are sometimes different from the entanglement witnesses for the density matrices 
under consideration. We stated and explored also a kind of "duality" where the 
partial transposition maps PPTflW to PPTflW, where W is another convex sub- 
set of LMM (the set of locally maximally mixed states), and the cases of bound 
entanglement detected in W are also cases for bound entanglement in W. 

Summarizing, we could present a detailed geometric structure of the subset 
of bipartite qudits under investigation. We think that this will help to find a 
good characterization of the whole state space and to investigate measures for 
entanglement for higher dimensional systems. 

In the outlook, we hope the paper advances our knowledge of these struc- 
tures, of the convex hull of higher-dimensional generalizations of the two-qubit 
Bell states. These two-qubit states have many applications in quantum infor- 
mation theory, and so their characterization for higher dimension is a desirable 
research goal. The higher dimensional generalizations also may have applications 
in quantum information, and display interesting geometrical features on their 
own. Since we show how to construct optimal witnesses explicitly and how to 
determine regions where there is bound entanglement, i.e. entanglement which 
cannot be distilled by local operation and classical communication (LOCC), these 
methods might be useful in quantum cryptography. Regarding possible applica- 
tions of our results and methods in a wider context, we note that the quasiclassical 
structure fits also exactly into the conditions needed for teleportation and dense 
coding, e.g. Refs. |BW92l IWOT] . 

Furthermore, according to the opinion of one of the referees, due to the Choi- 
Jamiolkowski isomorphism between a class of bipartite states and maps [ZB04J 
or channels [HHH98b] one may infer from the results of the present paper further 
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conclusions about the maps eg. (bistochastic) superpositive maps, which corre- 
spond to separable states: Using the isomorphism (for a recent exposition see 
e.g. [BQ6J or [L06] ) we see that the states of a d x d system satisfying one partial 
trace condition Tr^p = ojb represent stochastic maps (completely positive, trace 
preserving linear maps) while states satisfying simultaneously both conditions 
Tr aP = &b, Tr^p = uja, represent bistochastic maps. Hence the set LMM for 
bipartite systems is isomorphic to the set of bistochastic maps. See also [BDS01J 
for relations to remote state preparation. 

Last but not least, exploration of entanglement properties in still more detail, 
using the high symmetry of the chosen sets of states, seems to be the nearest and 
next goal. 
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